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i-C ' Abstract 

The commonly accepted belief that non- Abelian and Abelian mod- 
els are different because of the presence/absence of instantons and/or 
perturbative asymptotic freedom is analyzed from a historical perspec- 
^ ■ five. The presentation covers the major developments which brought 

about this dogma, as well as all the supportive evidence produced 
since. For a model possessing both asymptotic freedom and instantons 
it is shown rigorously that a disorder variable varies nonanalytically 
with the temperature. 



"'^This paper reflects the state of afi^airs as it was in 1991. In some areas there have been 
important developments since then. 
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1. Introduction 



The generally accepted belief among condensed matter and particle physi- 
cists is that there is a dramatic difference between Abelian nonlinear a- 
models in two dimensions {2D) and gauge theories in AD and their non- 
Abelian counterparts. This opinion, which did not exist prior to roughly 
1973, has its origin in two important observations: 

i) the importance of topological properties (Kosterlitz and Thouless |]I|]) 

ii) the discovery of perturbative asymptotic freedom (Gross and Wilczek 
§ and Politzer [|). 

The difference between the Abelian and the non-Abelian models is sup- 
posed to stem from the presence of instantons and/or asymptotic freedom in 
the latter, but not in the former. In the previous sentence we used 'and/or' 
deliberately, to emphasize that one property does not entail the other, hence 
one may legitimately ask, which property, if any, is responsible for this ac- 
cepted difference. 

We are not aware of any paper addressing this issue. Instead many papers 
report evidence in favor of what we shall call the orthodoxy, coming from all 
sorts of sources. For instance, it is argued that QCD4 must be the correct 
theory of strong interactions because there is general agreement between its 
predictions and experiment (for a recent example see 0]). In reality the 
confrontation of QCD4 with experiment is nowhere near that of say QED4, 
for two reasons: 

i) Any realistic computation must address the issue of constructing physi- 
cal states (confinement, etc), which even QCD believers admit is beyond their 
present technical abilities. Thus the best QCD predictions involve, besides 
some well defined perturbative computation, some ad hoc phenomenological 
assumptions. Hence the situation is not at all like computing the electron 
g — 2 factor in QED. 

ii) Even with all these uncertainties, the overall agreement between theory 
and experiment is at the 10% level, rather than 10~*, as in QED4. Under 
these circumstances, to employ nature as an analog computer for learning 
the true properties of QCD4 seems rather cavalier. 

Of course nowadays the computer offers the opportunity to investigate 
numerically the properties of many models of interest. Numerous such studies 
have been performed and the general consensus is that they do corroborate 
theoretical expectations. Finally many people have tried to produce indirect 
evidence in favor of the orthodoxy by appealing to so-called exact solutions. 
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the expansion, etc.. 

In this paper we offer a critical analysis of all the evidence supporting 
the idea that Abelian and non-Abelian models have fundamentally different 
properties. Without advocating the opposite point of view, we recall rigor- 
ously established facts and thus attempt to put all the supportive evidence 
in its proper perspective. After giving a historical overview of the subject we 
begin by recalling what is known about perturbation theory and semiclas- 
sical (instanton) approximations. We then review the more recent findings 
regarded as supportive of the orthodoxy, including Monte Carlo studies. Our 
general conclusion is that the original belief that the presence of instantons 
and/or asymptotic freedom explains the difference between Abelian and non- 
Abelian models is unfounded. The same conclusion applies to the more recent 
circumstantial evidence. Finally we argue that numerical studies are at best 
inconclusive - one may even claim that they indicate the opposite. Also in 
this paper we give an outline of a rigorous proof that a certain modifica- 
tion of the nonlinear 0{N) models in 2D, which retains asymptotic freedom 
and instantons, has a nonanalytic change in its behaviour with the inverse 
temperature f5 (a similar modification and result were proved by Mack and 
Petkova for gauge theories 0). 



2. Historical Background 

In 1966 Mermin and Wagner introduced their celebrated theorem 
about the absence of symmetry breaking in 2D models enjoying invariance 
under a continuous, compact group. Initially it was thought that that result 
implied that 2D nonlinear 0{N) cx-models cannot undergo phase transitions, 
since the theorem says that they could not exhibit long range order (l.r.o.). 
Shortly thereafter, however, Stanley and Kaplan examined the high tem- 
perature expansion of the magnetic susceptibility in such models. Although 
one cannot rigorously compute the radius of convergence of a series by know- 
ing only a finite number of terms, the numbers suggested very strongly that 
the magnetic susceptibility diverged at finite /?. To reconcile such a property 
with the absence of l.r.o., the term algebraic order was coined and initially 
it was believed that all 0{N) models with N > 1 possess such a phase at 
sufficiently low temperatures. 

This paradigm survived until 1973, when the seminal paper of Kosterlitz 
and Thouless appeared They introduced the notion of topological order; 



3 



specifically they argued that at low temperature, the typical configurations of 
the system will be low lying excitations of the configuration of lowest energy. 
For 0(2) they will consist of bound vortex-antivortex pairs, since the energy 
of just one vortex diverges as InL (L-linear size of the lattice); for 0(3) one 
can form instantons, whose energy is 0{L°). (We use the term instanton, 
which was invented only later, in 1976; the properties of the instanton con- 
figuration were described and employed by Kosterlitz and Thouless in Sect. 6 
of their 1973 paper.) Since the entropy of these topological defects is O(lnL) 
(the position of their center), in the 0(2) model, at low temperature, it 
will be overwhelmed by the energy and, Kosterlitz and Thouless argued, the 
model will enjoy topological order. On the contrary, in the 0(3) model the 
entropy will win and thus correlations will exhibit exponential decay even at 
low temperature. 

The proposal of Kosterlitz and Thouless pertained only to spin systems 
and was not readily applicable to particle physics. However the same year 
particle physicists became very excited by their own discovery, namely that 
of asymptotic freedom by Gross and Wilczek and Politzer . The proof 
that in perturbation theory non-Abelian gauge theories could enjoy this prop- 
erty promised a resolution to the long sought field theoretic explanation of 
strong interactions and even the grand unification of weak, electromagnetic 
and strong interactions (GUTS). Moreover it was claimed that asymptotic 
freedom suggests naturally why quarks are confined: if the interactions are 
getting weaker at short distances, they must be getting stronger as one tries 
to separate a quark from an antiquark, hence asymptotically free theories 
should also be confining p. 

For approximately two years the developments in condensed matter and 
particle physics remained largely separated. However in 1975-1976 in papers 
destined to become classics, Polyakov [H, Belavin et al |TD[ and Brezin and 



Zinn- Justin |TT|] bridged that gap. Namely Polyakov and Brezin and Zinn- 
Justin proved that the 0{N), N > 3, 2D nonlinear cr-models are asymp- 
totically free in perturbation theory, while Belavin, Polyakov, Schwartz and 
Tyupkin |]T^ showed that Yang-Mills theories in AD also possessed instan- 
tons - just as the 0(3) spin model in 2D. Through these realizations a new 
paradigm was born: there exists a fundamental difference between Abelian 
and non-Abelian models, stemming from their different topological properties 
and/or existence or absence of asymptotic freedom. 

This fifteen year old dogma remains unproven, yet it is widely believed. 
We will analyze its merits in the next section, but first we would like to finish 
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our historical recollection. During the late 1970s several attempts were made 
to develop a semiclassical approximation based upon the instanton idea in 
both 2D (T-models and AD gauge theories. In spite of the technical brilliance 
of these papers, Patrascioiu [|13| argued that the original computations by 
t'Hooft 0] (for gauge theories) and Berg and Liischer [|l5l and Fateev, Frolov 
and Schwartz |jl6[ (for a-models) were incorrect and that uncontrollable in- 
frared divergences plagued the semiclassical approximation. This fact was 
proven by Patrascioiu and Rouet for both gauge theories [jl^ and cr-models 
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These mathematical facts, established in 1981, remain largely ignored by 
the community, which still talks about the resolution of the U{1) problem by 
instanton effects or about the strong CP- violation problem caused by instan- 
tons. In reality the U{1) problem does not exist because the corresponding 
axial current has the famous Adler-Bardeen anomaly and nobody has ever 
shown that there exists some conserved, gauge invariant axial U{1) current 
in the physical sector. As for the 'strong CP problem', no such problem 
existed prior to 1976 and its origin is intimately connected with the belief 
that instanton effects force QCD4 to have many 6'-vacua; since all but one 
such vacuum would be CP violating, it was concluded that without 'fine 
tuning', the theory would not exhibit CP invariance. The axion was intro- 
duced in 1977 precisely to eliminate this need for fine tuning. Extensive 
experimental searches have repeatedly failed to find this elusive particle. Is 
this an embarassement for QCD^ ? Not in our opinion. The original instan- 
ton motivation was discredited by the infrared divergences found. Moreover 
the better understood QED2, which also has topolgically nontrivial fields, 
revealed that while '^-vacua' can be defined by choosing suitable boundary 
conditions, there is no necessity to have 6 different from [^, 

A rigorous result supportive of the orthodoxy was obtained the same year. 
Frohlich and Spencer |^ proved rigorously that Abelian models do undergo 
the phase transition predicted by Kosterlitz and Thouless. As far as direct 
evidence goes, this can be considered the only positive result ever found and 
everything else seemed to cast doubt rather than support the orthodoxy. For 
instance after a real tour-de-force Bricmont, Fontaine, Lebowitz, Lieb and 
Spencer proved that in Abelian models ordinary perturbation theory 
does provide the correct asymptotic expansion in powers oi 1/(5 for (3 ^ 00 
(the proof is for an infinite lattice model, so it does not concern constructing a 
continuum limit but only a thermodynamic limit). In spite of all their efforts 
(private communication to A. P. by J. Lebowitz), they could not extend the 
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proof to non-Abelian models. Of course the reader could consider this fact 
merely a measure of the prowess of mathematical physicists, however there 
are some exactly soluble cases, which we will discuss in the next section 
and they are not very reassuring - perturbation theory produces (infrared) 
finite answers, which are correct in Abelian models and false (at 0(l//9^)) in 
non-Abelian ones. 

Another potential embarassment was pointed out by Solomon |Q. By 
analyzing the high temperature expansion, he realized that if in the 0{N) 
model one changes the action from Si ■ Sj to {Si ■ SjY, the series seems to 
predict a singularity at a finite positive (3 for all < oo. This would be 
strange because this action leads also to asymptotic freedom for A^ > 3. 
Two numerical studies triggered by Solomon's work |]2S|, ^| observed a 



dramatic increase in the correlation length and magnetic susceptibility at 
some finite p. Although these two studies produced practically identical 
data, in a classic example that beauty is in the eye of the beholder, Fukugita 
et al concluded that a transition to a massless phase was occurring while 



Sinclair ||2^ claimed that it was just a cross-over regime. 

While one may easily dismiss numerical results, a rigorous result prov- 
ing that asymptotic freedom could not be the source of the mass gap was 



obtained by Richard ||2^. Let us state it for simplicity for 0(3): consider a 
modification of the measure such that sin 6 > e where 6 is the angle specifying 
the latitude. This modification destroys 0(3) invariance and instantons, but 
not asymptotic freedom. By using Ginibre's inequalities Richard proved that 
provided e > c/ ^/J3 for some suitably chosen c, for P sufficiently large the 
correlations of S^ and Sy must decay algebraically. Inspired by this result, 
Patrascioiu questioned the accepted relationship between the perturba- 
tive Callan-Symanzik /5-function and the nature of the spectrum, pointing 
out several counterexamples. 

As we stated in the introduction, the community has largely ignored all 
evidence shedding doubt upon the orthodoxy. Instead supportive evidence 
has been sought in many places. Since it is commonly argued that such 
supportive evidence has been found and that it is only a technicality that 
nobody has managed to prove rigorously that non-Abelian models are indeed 
fundamentally different from their Abelian counterparts, in the remainder 
of this paper we will analyze these findings in detail. We will also adapt 
the Mack-Petkova modification of gauge theories to 2D 0{N) nonlinear a- 
models. For these models we will prove rigorously that although they are 
asymptotically free and possess instantons a certain disorder variable changes 
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nonanalytically with (3. Were it not for the fact that disorder variables are 
nonlocal observables in terms of the original spins, this would constitute a 
rigorous proof for the existence of a phase transition. 



3. What is the Evidence for the Orthodoxy ? 

a) Asymptotic freedom 

Asymptotic freedom was shown only in perturbation theory which 



is an expansion in small deviations from an ordered state. As stressed by 
Patrascioiu and Richard the Mermin- Wagner theorem guarantees that 



such a state does not exist m. D <2. Since in spite of this difficulty it was 
possible to prove rigorously pSf that in Abelian models perturbation theory 
does provide the correct asymptotic expansion, one could take the optimistic 
point of view that perturbation theory does work also in non-Abelian mod- 
els and only technical difficulties are preventing a proof. However in ID, 
where one knows the exact answer, one sees by explicit computation that 
perturbation theory fails for non-Abelian models, while giving correct an- 
swers for Abelian ones. (In fact since in ID there is a nonvanishing mass gap 
for any finite /3, the infinite volume limit of the Green's functions is unique; 
yet in non-Abelian models, perturbative answers depend upon the boundary 
conditions used.) 

In gauge theories we expect similar difficulties with perturbative predic- 



tions. Indeed Patrascioiu pointed out |30] that even after the gauge freedom 



has been completely eliminated, on an infinite lattice the fiuctuations do not 
go to zero as /3 — ^> oo, as they do on a finite lattice. We would like to em- 
phasize that this result, proven in the complete axial gauge, is much stronger 



than Elitzur's theorem on the impossibility to break spontaneously a 
gauge symmetry and in our opinion strongly suggests that perturbation the- 
ory is incorrect in non-Abelian gauge theories; this can be verified in 2D 
(using periodic boundary conditions). 



h) Instanton computations 

We cannot delve here into the complicated instanton computations and 
refer the reader to the original papers quoted before (see also [^]). We will 
only try to explain from a heuristic point of view why there is a problem. It 
is well known that a semiclassical approximation involves two steps: 
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i) find a classical solution and 

ii) calculate Gaussian fiuctuations around it 

The second step amounts to the calculation of the determinant of an operator 



which in the simplest case of one instanton takes the form |jT^ 



O = -A^\-A + V{x)) = 1- A-^V{x) (1) 

Here V{x) is the operator of multipication by some instanton induced 
potential and the inverse of the Laplacian appears to ensure proper nor- 
malization of the functional integral. The determinant has an ultraviolet 
divergence even in a finite volume which can be cancelled by local countert- 
erms (see for instance |53|, M), giving rise to a 'renormalized determinant'. 



But here we are concerned with something else, namely an infrared diver- 
gence: since an instanton configuration has nontrivial topology, V{x) behaves 
cLS I X I cLS I X I ^ CX^ and thus it is an infrared singular perturbation of the 
Laplacian. Just as in nonrelativistic quantum mechanics Levinson's theorem 
fails for long range potentials, the renormalized determinant of such an oper- 
ator fails to exist in the infinite volume limit. In technical terms |Q, neither 
O nor any finite power of it are trace class (actually not even compact) and 
therefore if one computes the (renormalized) determinant for a sphere or a 
ball of radius R, the limit R —>■ oo does not exist since the determinant 
contains a term diverging like log R, as found explicitely by Patrascioiu and 
Rouet [00. 



c) 1/N expansion 

Over the years, many authors have used results obtained in the ex- 
pansion as evidence that perturbation theory does provide the correct asymp- 
totic expansion (see for instance and references given there). Let us use 
the usual scaling and write 



(3 = NP (2) 

In his well known paper Kupianen proved that the expansion is an 
asymptotic expansion at fixed f3. His error estimates are such that they do 
not allow to interchange the limits /5 — oo and — > oo, in particular not for 
long range quantities such as the correlation length ^. In fact the numerical 



data produced by Wolff ||3^, ^ ^ show very clearly that with increasing (3 
one has to go to higher and higher to achieve a certain degree of closeness 



8 



between the correlation lengths of the 0{N) model and the 0(oo) (spherical) 
model at a given /?. 

So there is no conflict between the successes of the expansion and 
the possibility that the 0{N) model undergoes a phase transition at some 
Pkt{N). What does follow however from Kupiainen's work (see the intro- 
duction of his paper) is that if (3kt{.N) < oo, then [3kt{.N) has to grow at 
least like a power of In for ^ oo. 



d) High Temperature Expansions 

Butera, Comi and Marchesini |5D, 5T| computed high temperature series 
(see also |^) up to order and Pade approximations to gain insight into 
the possible singularity structure of the 0{N) models in the complex /5-plane. 
They conclude that the 0(2) model has the nearest singularity on the positive 
real axis (and identify this singularity with the KT transition), whereas for 

> 3 they find that the closest singularities are off the real /3-axis. This 
is taken as supporting the absence of a phase transition in those models. It 
should be remarked, however, that their results do not provide any evidence 
against singularities on the real axis that are further from the origin than the 
complex conjugate pair they find. Pade approximants are notoriously unable 
to find such singularities which are 'shielded' by closer ones. 

Bonnier and Hontebeyrie used Pade resummation in a conformally 
mapped variable that is based on assuming the asymptotic scaling predicted 
by the perturbative /3-function. They report good agreement with Monte- 
Carlo data, but unfortunately the data they are using are very old ones of 
poor quality. The agreement deteriorates markedly if one is using the better 
data now available p8| . 

There is another point that should be noted: For N = 2,3,4 the high 
temperature series for the susceptibility has only positive coefficients up to 
the order to which it has been computed [|^. But a power series with positive 
coefficients has its nearest singularity on the positive real axis; so unless one 
believes that the coefficients computed so far somehow know already that 
the higher ones will eventually change sign, they cannot credibly predict an 
imaginary part for the closest singularity. 

So it seems fair to say that high temperature expansions, despite their by 
now remarkable length, have not produced any conclusive evidence in favor 
of the orthodoxy. 
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e) 'Exact solutions' 

Zamolodchikov and Zamolodchikov obtained an 'exact 5'-matrix' for 
the continuum 0{N) models under a number of assumptions: First of all 
they assumed that the model describes an 0{N) vector multiplet of massive 
Bose particles; furthermore they made the usual assumptions about unitar- 
ity, analyticity and crossing symmetry, and finally the less standard ones 
of absence of particle creation, minimal singularity structure and most im- 
portantly, factorization. (Absence of particle creation and factorization are 
supposed to follow from the existence of infinitely many conservation laws 
p5| , whose existence, however, in turn depends on some assumptions 
and cannot be proven without a construction of the continuum limit of the 
model). It is clear that their construction, remarkable as it is, cannot help 
answer the question of the existence of a massless phase of the lattice models, 
since it assumes a mass gap from the beginning. 



Using similar assumptions, Karowski and Weisz derived an 'exact 
current form factor', that is the matrix element of the current operator be- 
tween the vacuum and 2-particle states. The same remarks as above apply 
to this construction. 

Polyakov and Wiegmann [^] produced an 'exact solution' of the 2D 0(4) 
model; Wiegmann |^ extended the method to the 0(3) model. The first 
step in this approach is to map the nonlinear a-model into a model with 
4-fermion interaction which then is to be solved via the Bethe ansatz. The 
problem is that to establish this equivalence, one has to use an identity for a 
Gaussian integral over a certain gauge field and ignore the fact that actually 
the gauge fields vary over the compact space SU{2) = S3 and the integration 
has to be done with the Haar measure, not the flat (Lebesgue) measure. Thus 
the steps required in going from eq.(l) to eq.(2) of cannot be justified on 
a lattice and are valid only if one imagines taking a naive continuum limit. To 
quote from a recent paper ||50| it is therefore not clear whether the solution 
'in addition to being exact, is also correct', and if so, for what model. 

Hasenfratz, Maggiore and Niedermayer compared the Polyakov- Wiegmann 
solution that depends on a mass parameter with the perturbation expansion 
and derived a formula for the mass gap of the 0(3) and 0(4) models; later 
Hasenfratz and Niedermayer carried out a similar calculation for general 
0{N) starting from the Zamolodchikov^ ^-matrix and comparing with 
perturbation theory. This formula has the property of giving in the limit 
N 00 the correct asymptotic behavior of the mass gap of the spherical 
model for (3 00. So it should come as no surprise that numerical tests 
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(see f) below) showed good agreement with the formula for large but 
considerable deviations for smaller A^. The derivation itself involves some as- 
sumptions that may be questioned, such as the description of a very dense gas 
of particles in terms of a 2-particle S-matrix. In addition, as already noted, 
the S'-matrix used as an input can only be considered a clever guess, consid- 
ering the many assumptions that are needed to obtain it, and its connection 
with any 0{N) lattice model is not in the least transparent. So whatever its 
merits, this work also does not contribute to answer the question of the mass 
gap, since its existence has to be assumed from the start. 



f) Numerical results 

Numerous papers have appeared reporting Monte Carlo investigations 
that are claimed to support the absence of a phase transition in the standard 



nearest neighbor action (s.n.n.a.) non-Abelian ferromagnets [|53| , |5; 



57|. We believe that to a large extent these claims were motivated by the 



authors' expectations and that in fact an objective analysis of the numerical 
situation suggests rather the contrary. Namely, there is universal agreement 
that in almost all the standard action models there is a 'cross-over' region, 
where the magnetic susceptibility and the correlation length increase faster 
than the asymptotic freedom predictions. This 'cross-over' region is supposed 
to reflect the existence of a line of first order transitions, terminating at a 
critical point in the 'mixed action models' parametrized by (/5i,/?2); where 

H{(3uf32) = - j:{f3iS{x) ■ S{y) + (32{S{x) ■ S{y)f} (3) 

We think that this scenario is highly implausible; indeed although not 
rigorously proven, one would expect that in a ferromagnet susceptibility and 
correlation length are nondecreasing functions of f3i and /32- Hence if they 
diverge at the point {Pi^kt, P2,kt), they must continue to do so in the whole 
region Pi > Pi^kt and P2 > P2,kt- But there are no phase boundaries 
separating this region from the line = 0. Hence the critical region must 
touch the line /32 = 0. Therefore we believe that there is no good explanation 
for the repeated occurrence of the so called 'cross-over' region observed in 
both 2D 0{N) models and AD gauge theories. On the other hand, it could 
very well be that this region is not a cross-over, but rather the neighborhood 
of a critical point, which the orthodoxy claims should not exist. 
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We would like to clear another erroneous belief expressed in many numer- 
ical studies. Several authors Q have advocated going past the 'cross- 
over' region by employing the Monte Carlo renormalization group (MCRG) 
or some finite size scaling curves. Common to all such approaches is the 
belief that one can take measurements on small lattices and learn about the 
infinite volume behavior. As we explained in (a) above, perturbation theory 
is suspect precisely because an ordered state does not exist on an infinite lat- 
tice. On the other hand, given the size of the lattice, one can always choose 
a (3 sufficiently large so that all these lattice Green's functions agree with 
their perturbative values to any degree of accuracy (since perturbation the- 
ory is clearly asymptotic in a finite volume). We have checked that already 
for (3 in the crossover region the finite volume susceptibility computed via 
Monte Carlo agrees with the perturbative formula two-loop formula given 
by Hasenfratz within a few percent, provided the size of the lattice is 
less than the infinite volume correlation length, and the agreement is rapidly 
improving with increasing (3. So there is really no insight to be gained by 
running Monte Carlo simulations in this regime. The real dilemma is, do the 
limits L — > cxo and 13 ^ oo commute? Techniques such as MCRG not only 
cannot answer this question but, by insisting on working on small lattices, 
are bound to reproduce all perturbative predictions. 

Liischer and Wolff ^D|] studied numerically the current form factor and the 
2-particle scattering phases of the 0(3) model and claimed agreement with 
the results of Zamolodchikov^ and Karowski and Weisz [^], respectively. 
They worked at values of (3 where the model is clearly in its massive phase and 
the correlation length is between 6.9 and 13.6 lattice units. The form factor 
is fixed by a normalization condition at zero momentum; for increasing mo- 
menta they find increasing deviations from the predicted values. According 
to them the differences can be understood as lattice corrections (a derivation 
of those corrections is not given). They also find that the 2-particle scat- 
tering phases roughly agree with the predictions of for energies that are 
small compared to the mass. Discrepancies beyond the numerical accuracy 
are again blamed on 'lattice artefacts' (because these discrepancies were in- 
creasing with the energy and were smaller at the smaller value of the two 
values of the mass investigated). It should be remarked that the numerical 
determination of the scattering phases depends on quite a lot of additional 
theoretical input and also on an extraneous parameter (called to). But it is 
most important to realize that by its very nature such a test cannot answer 
the question of the existence of a mass gap, and due to the limitations of nu- 
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merics it could also not verify the interesting feature of the Zamolodchikov^ 
S'-matrix that the high energy limit of the scattering phases is zero, which 
has been interpreted as a manifestation of asymptotic freedom. 

Wolff ||3^ performed a numerical study to test the validity of the mass 
formula of Hasenfratz and Niedermayer |^ mentioned above for = 3, 4, 8, 
of course in a regime where there is undoubtedly a mass gap. It turns out the 
the agreement is not very good for = 3 (between 25% and 33% discrepancy 
in the region beween (3 = lA and 1.9), but seems to improve with increasing 
A^. In our opinion this study has no bearing on the question of the existence 
of a massless phase, since the data were taken in the same range of /3 = 
(3/N for the different values of A^. Of course for A^ — oo at fixed and 
sufficiently large (3 the agreement has to improve because, as remarked above, 
the Hasenfratz-Niedermayer formula is correct for the spherical model in the 
limit /3 — cxD, and at fixed /? the mass gap of the O(A^) model seems to 
converge to the spherical model mass (Kupiainen ||36| proved that the mass 
gap is bounded below by a quantity converging to the spherical model mass). 
But we know anyway from Kupiainen's work |^ that (3kt has to increase at 
least logarithmically with A^. So only a study of the mass gap in such a regime 
could possibly give any information about the existence or nonexistence of a 
massless phase; this is beyond the present numerical possibilities even with 
the new cluster algorithms. 



g) The role of topology 

The original Kosterlitz-Thouless scenario was that in the 0(2) model a 
phase transition must occur, refiecting the loss of topological order as the 
temperature is increased. Their conclusion was derived from the energy- 
entropy arguments mentioned earlier: on a lattice of size L a vortex has an 
energy of order log L. Its entropy, measuring essentially the location of the 
center is also O(logL). Hence if [3 is too large, vortices are bound, while 
at high temperature they unbind, triggering the phase transition. These 
considerations suggest also a basic difference between 0(2) and 0(3); in the 
latter smooth configurations - instantons - have energies 0(L°), hence they 
act like point defects and disorder the system at arbitrarily low temperatures. 

To understand better the role of topology, we found it useful to consider 
a modification of the 0(2) model dubbed 'cut' model in ||6T[, 'constrained 
model' in ||62|: the Gibbs factor retains its form only for \S{x) — S{y)\ < e, 
while it is replaced by otherwise. Thus we are forbidding large angular de- 
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viations between neighboring spins. Since this modification is ferromagnetic, 
we would expect (3kt to decrease as e is decreased from its s.n.n.a. value 
of 2. This is indeed what we observed numerically. The surprising result 
though was that we found that for about e = 1.57, Pkt = on a square 
lattice; for this value of e vortices are still allowed and at /3 = they cost no 
energy! This finding suggests that the original Kosterlitz-Thouless argument 
was too naive in its estimate of the entropy: at /3 = 0, the Gibbs factor is 
either 1 (for \S{x) — S{y)\ < e) or (otherwise). So on a square lattice, the 
Kosterhtz-Thouless argument would have suggested that Pkt = only for 
e < \/2, since that is the value for which vortices cease to exist. 

The low temperature phase of the 0(2) model must thus be characterized 
by the fact that vortices are sufficiently rare and spin waves dominate. It 
was suggested already by one of us in ||63| that in both Abelian and Non- 
abelian models at large (3 spin waves may dominate and defects may be 
suppressed. In one of us is giving more detailed arguments in favor of 
this scenario. According to this scenario the situation would be as follows: 
just as in > 3 the very good 'local order' present at large P manifests 
itself thermodynamically as l.r.o., in 2D its manifestation may be algebraic 
order. As the temperature is raised, defects - bonds where \S{x) — S{y)\ is 
large - become more abundant and at some point condensate, putting the 
system in a phase characterized by exponential decay of correlation functions. 
Artificially suppressing defects should have an effect in all dimensions and we 
predict that a suitably modified 0{N) model will exhibit l.r.o. at all /3 > 
in any D > 3. 



h) QCDi 

One may wonder what this discussion implies for QCD. It has generally 
been accepted that there is a close analogy between 2D spin models and AD 
gauge models. Most of our discussion applies there equally well, with the 
proviso that the evidence for the orthodoxy is weaker: the 1/A^ expansion 
does not have a firm base like Kupiainen's ||3^ results, there is no exactly 
soluble N = oo limit, there are no 'exact solutions' for the S'-matrix, the 
role of topology for disordering the system is even less clear and finally the 
numerical results are much more limited by the increased requirement of 
computing power. 

The existing numerical data certainly cannot rule out a deconfining phase 
transition and the famous 'dip' in the /5-function may in fact be an indication 
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of its presence, since it shows that quantities hke the string tension tend to 
zero faster than asymptotic scahng would predict, just as in the s.n.n.a. 0{N) 
models in the so-called cross-over region. The usual explanation, that the 
'dip' is a cross-over regime induced by the presence of a critical point in the 
Pfund — Padj plaue |^ , could again be ruled out by correlation inequalities, 
hence seems implausible. 

Finally, let us recall the following bizarre situation: it is known that all 
gauge theories undergo a deconfining transition at finite temperature ||66 . 
It is also agreed that U{1) and SU{N), N > 4 gauge theories undergo a 
transition at zero temperature. For the Abelian case U{1) that transition 
is known rigorously [^, 37] to be deconfining. For the non- Abelian cases 
SU{N), N > 4, the claim is that the transition is first order and not decon- 
fining |6^, 



73, 741 that it is difficult 
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It has been recognized 
to separate this first order 'bulk' transition from the finite temperature de- 
confining transition, because the strange thing is that the extrapolation of 
the curve describing the deconfining transition in the (3 — T plane (T is the 
temperature) seems to cross the axis T = just where the supposed first 
order transition occurs. In the orthodox outlook, these two transitions are 
supposed to have nothing to do with each other (in fact we know of no expla- 
nation for the first order transition at zero temperature, which for instance 
is not seen in SU{2) and SU{?>)). Then one could ask: if in SU{N), N > 4, 
these two transitions have nothing to do with each other, should it not be 
possible to change the action and pull the two transitions apart, perhaps 
even eliminating the first order transition altogether? Could it not be that 
the difficulty of separating the two transitions simply means that the first or- 
der transition at zero temperature is nothing but the deconfining transition? 



This possibility was discussed in [75 



4. Rigorous results for a modified 0(N) model 

To illustrate that the existence of instantons and perturbative asymptotic 
freedom do not rule out nonanalytic behavior we consider the following model 
which is inspired by some considerations for lattice gauge theories due to 
Mack and Petkova 0: It is defined by the following modification of the 
s.n.n.a. measure: for every plaquette p, we require 
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n sgn{S{x)-S{y)) = l (4) 

{xy)&dp 

where the product is over the bonds around the plaquette p. 

Remark. This modification should be unimportant at low temperatures. 



Indeed, using a result of Bricmont and Fontaine the probability to vio- 
late the constraint at a given plaquette in the s.n.n.a. 0{N) model is bounded 
by exp(— c/3) for some c > 0. Moreover the modification does not affect the 
existence of asymptotic freedom in perturbation theory nor the presence of 
instantons (defects costing an energy of order 1) in 0(3). Finally we remark 
that the 'cut' or 'constrained' models discussed in |6l|, |6^ (see also point h) 
of the previous section), in which the angular deviation between neighboring 
spins is limited, automatically satisfy the constraint (4). We note that the 
constraint should make the system more ferromagnetic so that one would 
expect 

{SiO)-S{x))m> {S{0)-Six))s.n.n.a. (5) 

where {.)m denotes the expectation in the modified model. For this model 
one can introduce the disorder parameter 

L-l 

{D{x))rn = lim (TT exp{-2Psgn{S{xk) ■ S{xk+i))))m (6) 

where the product is over a path (ordered set of bonds) from x to x^. Because 
of the constraint (4) {D{x))m is path independent. The following theorem 
can be proven along the lines of Mack and Petkova ||^: 
Theorem: 

a) There exists a /?i G (0, oo), such that for P < Pi 

{Dix))m > 

b) There exists a P2 ^ {Pi, 00) such that for P > P2 

{D{X)):^ = 

Remark: The only reason why one could doubt that this theorem implies the 
existence of a phase transition in these 0{N) models is the non-local nature 
of the disorder variable D{x). But the general experience gathered in the 
study of numerous models is that a disorder operator of the kind used here 
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signals correctly the occurrence of a phase transition and can even be used 
to determine exactly its location. 

Proof: We do not give a detailed proof, because the adaptation of the 
proof given by Mack and Petkova for gauge theories is straightforward. We 
will only present the strategy. 

Because of the constraint (4) one can define Ising spins {cTx} attached to 
the sites as follows: 

L-l 

c^x = n sgn{S{xk) ■ S{xk+i)) (7) 

A:=0 

where as before the poduct is over a path (of length L) from to x and by 
(4) there is no dependence on the path. The system can now be rewritten 
as a coupled system of the 0{N) spins S and the Ising spins a, described by 
the Hamiltonian 

H=J2 ^o:ayS{x) ■ S{y) (8) 

with the constraint that for all bonds {xy) S{x) ■ S{y) > 0. As in 
one carries out a duality transformation of the Ising variables, leaving the 
5* variables alone. One obtains again a certain Ising model with fluctuating 
coupling constants given by 

Pxy = -^lntanh(/?(^(x) ■ S{y)) (9) 
Then as in the GKS correlation inequalities imply that 

{D{x))m> \{D{x))is\ (10) 

where the right hand side denotes the expectation value of the standard disor- 
der operator in the Ising model at the same value of p. (10) simply expresses 
the fact that the coupling of the Ising variables in the model described by 
(8) is less ferromagnetic than in the standard Ising model at the same f3. So 
we have learned that for f3 > (3c,is = |ln(-\/2 + 1) {D{x))m = 0. In other 
words the first part of the theorem is proven for 

/5i>^ln(V2 + l) (11) 



17 



As is well known, the disorder parameter will become the two-point func- 
tion of the Ising spins after the duality transformation. The crucial point 
is that jSxy will typically be very small, in other words, the dual model will 
be in its high temperature phase, and its 2-point function will decay expo- 
nentially. The actual proof proceeds by cluster expanding the conditional 
expectation value of the dual Ising two point function for fixed couplings 
{Pxy}- One takes the remaining (annealed) expectation value over these cou- 
plings termwise in the cluster expansion. These expectations can be bounded 
above using the so-called chessboard bounds (after undoing the duahty trans- 
formation), since the measure still possesses reflection positivity. In this way 
one obtaines convergence of the cluster expansion and exponential decay of 
{D{x))m- For details of this proof we refer the reader to the paper of Mack 
and Petkova 0. 

We would like to remark that something more is true for this model: 
the Ising spins (7) of this model show a transition from a disordered high 
temperature phase to a phase with l.r.o.. The first part of this remark is a 
direct consequence of the GKS inequalities. To prove the second part one 
has to adapt Georgii's result on the percolation of the low energy bonds 
to this model; this result implies in particular that the bonds (xy) with 
sgn{S{x) ■ S{y)) > percolate. Using the definition (7) one sees that this 
implies the existence of a percolating cluster of + Ising spins, implying long 
range order. 
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